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Abstract 

We analyze the dynamical evolution of the resonant tunneling of an ensemble of 
electrons through a double barrier in the presence of the self-consistent potential created 
by the charge accumulation in the well. The intrinsic nonlinearity of the transmission 
process is shown to lead to oscillations of the stored charge and of the transmitted and re- 
flected fluxes. The dependence on the electrostatic feedback induced by the self-consistent 
potential and on the energy width of the incident distribution is discussed. 
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In recent years there has been renewed interest in the phenomenon of resonant 
tunnehng (RT) through double barriers. The unique capabilities of molecular beam epitaxy 
make it possible to investigate fundamental questions on RT through simple man-made 
potentials by controlling the barrier and well parameters (e.g. height, thickness or barrier 
phase area) down to the atomic scalet-*^!. 

In this paper we investigate the dynamics of RT of ballistic electrons in the presence 
of the potential created by the charge trapped within the barriers This problem is 
interesting not only from a technological point of view but also as a test of quantum 
mechanical non-equilibrium situations in which many particles are involved. 

The model we propose tries to describe the following situation. A bunch of elec- 
trons is created within a contact layer and launched towards embedded semiconductor 
layers forming a double barrier potential. The charge dynamically trapped by the reso- 
nance will produce a reaction field which modifies the time evolution of the system. An 
exact treatment of such a problem looks very complicated. We assume a decoupling be- 
tween the longitudinal (in the direction x of motion perpendicular to the double barrier) 
and transversal degrees of freedom. This is a common assumption in treating tunneling 
phenomena. It makes the problem one-dimensional and allows the following factorization 
of the wave function: 

*(xi,X2,...,xjv;t) = ^{x\,X2, ...,XN;t) i^{yi,zi,y2,Z2,...,yN,ZN;t) (1) 

where x = x,y,z are the coordinates of each electron and t the time. We remark that 
antisymmetry of ^ is established if we take, for example, ^» symmetric and fl antisymmetric 
in their arguments. The experimental set-up we have in mind puts ideally all the electrons 
in the same high-energy longitudinal state while the transversal degrees of freedom are 
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essentially decoupled I^l. Therefore the choice (1) with $ symmetric is the only possibility. 
In other words the transversal degrees of freedom assure the respect of Pauli principle. 

Finally we assume the electrons in the bunch at the initial time uncorrelated which 
corresponds to a choice of $ as a product state of single-particle states iIj{x,0). At this point 
theorem 5.7 of Ref. [3] guarantees in the mean field approximation (which is reasonable 
due to the large number of electrons involved) that the state $ remains a product state 
during its evolution and allows us to write the following self-consistent equation for 'ip{x,t): 



+ 



'^{x,t) (2) 



2m dx"^ 

The external potential, V{x), is assumed, as customary, to be a step function (no 
electric field is applied): 

V(x) = Vo[d{x - a)9(b - x) + 9{x - c)9{d - x)] (3) 

with a < b < c < d and where 9{x) is the Heaviside function. The nucleus W{t, t'; x,x') is 
modelled assuming that memory effects can be neglected, i.e. W oc S{t — t'). We represent 
the global repulsive feedback effect, induced by the charge localized in the well, by a shift 
of the bottom of the well to a higher energy, Vqit)^ proportional to the charge, i.e.: 



/ 



W{t, t'; X, x')\ijj{x', t')\^dt'dx' = aVo^^9{x - b)9{c - x) (4) 

Qo 



Q{t) is the charge localized in the well at time t and Qq is a normalization charge which 
depends on the shape of the initial state, assumed localized around xq: 

/■c i-xo + ic—b)/2 

Q{t)= / dx'\^{x',t)\^ Qo= / dx'\iP{x',0)\'' (5) 

Jb Jxo-{c-h)/2 

Qo introduces an artificial dependence of the Eq. (2) on the initial condition. We have used 
this parameterization to make the comparison of different numerical simulations easier. 



The parameter a in Eq. (4) can be varied to reproduce phenomenologically the response 
of the medium to the charge trapped in the well and the characteristics of the electron 
bunch, e.g. its areal density. 

The 1-particle state which is the initial condition in our mean field equation (Eq. 
2) has been chosen to be a gaussian shaped superposition of plane waves with mean mo- 
mentum fikQ-. 



0) = i— exp 



1 / ^ 2 

1 ( X — Xq 



2 V a 



+ ik^x 



(6) 



with energy spread (energy full width at half maximum of the square modulus of the 
Fourier transform of (6)) Fq = 2\/ln2 h^ko/ma. xq is chosen so that at the initial time no 
appreciable charge sits in the well, i.e. (5(0) = 0. 

The solution of the differential equation (2-5) with the initial condition (6) has 
been achieved by a numerical integration on a two-dimensional latticet^^. Assuming for the 
barrier and well widths the values b—a = d—c = 20 oq and c—b = 15 oq (ao — 0.529 A being 
the Bohr radius) and for the barrier height 0.3 eV and using for m the free electron mass, 
the potential V{x) exhibits a single resonance in the transmission coefficient at energy 
Eji ~ 0.15 eV, the shape of which is well approximated by a lorentzian of full width 
Fr ~ 5 meV. The choice of these parameters was a compromise between the requirement 
of standard technological values and that of reasonable computation times. The electron 
mass was set to its free value to avoid the complications of a space variable effective mass. 

The incoming state mean-energy has been chosen so as to satisfy the resonance 
condition h^kQ/2m = Er. Finally, the normalized charge in the well, Q(t)/Qo, has been 
plotted as a function of time. The results for different choices of the parameter a are shown 
in Figs. 1-4 in the case of states with energy-spread larger, of the order of and smaller 



than the resonance width. 

The Fig. 1 represents the scattering of a wave packet on the fixed double bar- 
rier (hnear Schrodinger equation). When the packet is energetically much wider than 
the resonance, the building up and the decay of the charge are asymmetric, the decreas- 
ing following!^] the law exp{—t/r) with r = fi/Tn (decay of a lorentzian-shaped quan- 
tum state). On the other hand, for a wave packet narrower than the resonance, the 
charge presents a symmetric behavior, like the law exp[—{{t — to)/''')^]) where r = a/vQ = 
2\/ln2 h/To, to ~ \{b + c)/2 — xo\/vo and vq — hko/m (free evolution of a gaussian-shaped 
quantum state for negligible time-spreading). The latter result is not surprising, since 
now the packet traverses the double barrier almost undistorted. In the case of a state 
energy-spread comparable with the resonance width, the evolution of the trapped charge 
interpolates between these two extreme behaviors. 

When the non-linear term is effective, i.e. a 7^ 0, the evolution of the trapped 
charge changes drastically and oscillations can appear. This phenomenon has been quali- 
tatively anticipated by Ricco and Azbel'^^ . Their reasoning was very simple. At the initial 
time no charge is present in the well, the electrons are moving towards the double barrier 
and the resonance condition is fulfilled. When some charge penetrates into the well, the 
modification of the potential destroys the resonance condition. As a consequence the quan- 
tity of trapped charge has a maximum followed by a decrease. The resonance condition 
tends to be restored and a new cycle begins. However, as it will appear in the following, the 
nonlinearity makes the interpretation of the phenomenon significantly more complicated. 
For example, the conclusion by Ricco and Azbel that the above effect should be maximal 
for monochromatic states, is not correct. 

A detailed analysis of Figs. 2-4 (that are a sample of our global numerical work) 
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suggests the following observations. Oscillations are present, for appropriate values of the 
strength of the non-linear term, a, only when the energy spread of the state is wider or 
comparable to the resonance width. No oscillations are seen for nearly monochromatic 
states. When a increases, the oscillations, if present, tend to increase in number but 
decrease in amplitude. 

To understand these results, let us first interpret the dependence of the intensity 
of the trapped charge as a function of the parameters a and a. We simplify the question 
by considering a time-average of the charge dynamically present inside the well. Since 
during the time evolution VQ{t) and Q{t) are related with each other by the Eq. (4), a 
similar relation has to hold between the relevant time averaged quantities denoted by Vq 
and Q. Let us suppose, now, that we have a time independent situation with the bottom 
of the well at level Vq. As can be shown by explicit calculations, the charge Q present in 
the well is a fraction, 7, of the asymptotically transmitted charge Qt^^^' 



where 'ip{k,0) is the Fourier transform of (5) and |tvQ(/c)P is the transmission coefficient 
of the depicted potential. The time-average of Eq. (4) can be combined with Eq. (7) to 
obtain a self consistent relationship for Vq (or Q): 



The two sides of this equation are plotted in Fig. 5 for different values of a and cr; their 
intersection points represent our estimate for the time-averaged normalized charge trapped 
in the well during the interaction of the packet with the double barrier. The factor 7 is 
fixed by imposing that for ck = the results of Fig. 1 are reproduced. As expected it is 




(7) 



Vq _ 7Qt{Vq) 



(8) 



oiVo Qo 
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of the order of unity. Fig. 5 predicts correctly the time-averaged obtained from Fig.s 2-4. 
We also notice that for very large a the charge tends to disappear due to the smallness of 
the transmitted amplitude. 

We then try to understand the oscillating behavior. Let us assume that this 
phenomenon is due to the competition of two processes: (a) the filling up of the well by 
the incoming wave packet and (b) the natural decay of the trapped charge. For the process 
(a) the time scale is of the order of h/Fo (this estimate appears more reasonable for states 
narrower in energy than the resonance width). For the process (b) a reasonable time scale 
is h/TvQ, where Fvg is the energy spread of the function to be integrated in Eq. (7) 
(spectral decomposition of the charge present in the well). Oscillations are then expected 
if a substantial crossover of Fy^ and Fq is realized for the Vq values reached during the 
time evolution. The analysis of the function |'0(/c, 0) Pl^Vg (^) P shows that Tvq rises from, 
approximatively, TqTji/ y^r^~+T^ at Vq = 0, to a maximum greater than Fq (position and 
amplitude of the maximum are roughly proportional to Tq/Tr) and, eventually, decreases 
to Fq. As a consequence when Fq ^ F_r, Fvq is very close to Fq and independent of 
Vq. No oscillations are possible in this case for any value of a. On the other hand, when 
To > Fit, crosses Fq at a some Vq; oscillations are then realized for a sufficiently high 
value of a. This critical value of a increases with the ratio To/Tn. These predictions agree 
quantitatively with the results of the simulations reported above. 

The effects we have discovered exhibit a considerable stability under modifications 
of the model for example in the direction of a more realistic choice of the self-consistent po- 
tential. A deeper analysis of the problem, together with an investigation on the possibility 
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of an experimental study, is matter for a forthcoming publication. 
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Figure Captions 



Fig. 1 Time development of the normalized charge trapped in the well in the case of 
linear Schrodinger equation, i.e. a — 0, for states with energy spread much wider 
(cr = 110 ao and Tq = 43.2 meV), of the same order of (cr = 825 ao and Tq = 
5.8 meV) and much smaller (cr = 5775 ao and Tq — 0.8 meV) than the resonance 
width {Tji ~ 5 meV). An atomic unit of time corresponds to 4.83 lO"-*^^ seconds 
and the Bohr radius ao is 0.529 A. 

Fig. 2 As Fig. 1 but in the case of effective non linearity in the Schrodinger equation 
with a = 0.1. 

Fig. 3 As Fig. 2 but with a = 1. The case with a = 5775 ao is not shown. 
Fig. 4 As Fig. 3 but with a = 10. 

Fig. 5 Self-consistent estimate of the time-averaged normalized charge trapped in the well 
for different values of the parameters a and a. 
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